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1 INTRODUCTION 

The GAIA mission will offer the wonderful opportunity to establish galactic dynamical models 
based on the data of a large number of stars. Indeed, complete positional and kinematical data 
over a large volume in the Galaxy are needed to construct well constrained models: we do not 
have such data for the moment, but this is precisely what the GAIA satellite will provide. 

We present here new tools to establish axisymmetric equilibrium models of the Milky Way, 
based on GAIA-like data. We know that our Galaxy is a barred one and so, the axisymmetric 
hypothesis is not correct, but axisymmetric models are a prerequisite for perturbation theory 
(for example the theory of spiral density waves) and thus they are still very useful. 

The models we wish to establish are pairs ( V , F) where V is the gravitational potential 
generated by the whole mass distribution including the dark matter, and F is the distribution 
function in phase space for late-type tracer stars in the galactic disk. For an equilibrium 
model (stationary potential and stationary distribution function solution of the collisionless 
Boltzmann equation), we know that the distribution function in phase space depends only on 
the integrals of the motion. During the last decade, many studies (e.g. [0) have shown that 
the distribution function of tracer stars in the Milky Way has to depend on three isolating 
integrals of the motion, especially if information on the vertical motion of the stars is present 
as it will be with GAIA. In order to have an analytic third integral, in addition to the binding 
energy and the vertical component of the angular momentum, we use Stackel potentials ||. 

2 THE POTENTIALS 

In this contribution, we continue the work of Batsleer & Dejonghe [||, who presented a set 
of simple axisymmetric Stackel potentials with two mass components (halo and disk) and a 
flat rotation curve. The existence of a thick disk as a separate stellar component is now well 
documented (e.g. Q): so we have generalized the Batsleer’s potentials by adding a thick disk 
to them. Our new Stackel potentials are described by five parameters. They turn out to have 
an effective bulge because the mass density grows faster than an exponential near the center in 
the plane. The other motivation to continue Batsleer’s work is that, in recent years, Hipparcos 
data have enabled an accurate determination of some fundamental galactic parameters in the 
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Fig. 1. Left panel: The mass density in the plane. Right panel: For a fixed value of the axis ratio of the thin 
disk and of the halo, the figure displays all the values of the axis ratio of the thick disk and of the contributions 
of the disks to the total mass that are consistent with Hipparcos observables. 


solar neighbourhood: the mass densityp© ]|,jG} and the Oort constants A and B jTj]. We have 
looked for the Stackel parameters that are consistent with the above observables: as can be 
seen on Fig.l, many different combinations of the parameters are found to be consistent with 
those observables. 

3 THE MODIFIED FRICKE COMPONENTS AND THE MODELING 

In the future we will test the ability of these potentials to be combined with a distribution 
function, in order to reproduce kinematical data like those that will be provided by GAIA. This 
distribution function will be expressed as a linear combination of basis functions depending 
on a few parameters. We have defined new component distribution functions (modified Fricke 
components) with 9 parameters, that depend on three integrals of the motion , and that can 
represent realistic stellar disks when a judicious linear combination of them is chosen in a 
realistic galactic potential. These components have a finite extent in the ^-direction and have 
realistic scale lengths. More details can be found in ||. 

To construct our model, the method will be iterative: we will choose a potential, find a 
linear combination of modified Fricke components that fits the data by using the quadratic 
programming technique described in and then modify the potential in the light of the 
quality of that fit. This will provide some new constraints on the mass distribution in the 
Galaxy and some information on the dynamical state of the different late-type stars (since the 
Galaxy is a system with a long memory, this could provide information on the history of the 
Milky Way). 
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